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Supporting Information 

In the following we use square brackets to refer to equations from the main text of the 
paper. 


I. STABILITY DIAGRAM FOR DOUBLE DROPS. 


In Fig. 1 we plot contour lines on the stability diagram corresponding to the condition 
cos 9j = const = c for various values of the constant c in the case of a double droplet (eq. 
[3]). One can clearly distinguish the limiting cases mentioned in the main text: 

(i) 7 i —* 7 ij and 'jj/'jij —■► 0. In terms of c eq. [5] can be rewritten as 


ijhij = -{aha ~ !)■ 


( 1 ) 


In this limit different c correspond to straight lines on the stability diagram with 
different slopes but all converging to a single point (77/7^ = 1, 7 j/ 7 *j = 0)- For 
j = A, B see Figs, la, lb, respectively. 


(ii) 7 ^ —> 7 i and (7* — 7y)/7p —»► 0. Eq. [6] can be rewritten in terms of c as 


Jj/lij Ji/lij C- 


( 2 ) 


Accordingly, in this limit different c correspond to straight parallel lines with the same 
slope equal to 1 but with different vertical shifts c. 


Outside the limiting linear regions there is a transition region (approximately for n /i/'Yij E 
[0.2,4]) in which the lines cos 0j = const are not straight. 


II. CALCULATION OF THE CURVATURE OF THE INTERNAL INTERFACE 
OF PERFECT JANUS DROPS. 

In the case of three-dimensional double Janus drops the line k*(0B) separating the regions 
of positive and negative curvatures iU 1 in the (9b, k) plane (solid line in Fig. 2 a) is given 
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Figure 1: (a) Lines of constant cos 9a = c (eq. [3]) for c = 1, 0.75, 0.5, 0.25, 0, —0.25, —0.5, —0.75, —1 
(from top to bottom), (b) Lines of constant cos 9b = c (the same values as in (a) but from bottom 
to top). Note that in both cases in the limit of large 7^4 and 7 b the lines are parallel and have 
equal spacing due to the constant increment Ac = 0.25 between the neighboring lines. 


by the following equation (analogical to eq. [10]): 

Mb) 


k*(9 B ) = 


1 — /o(0b) 


( 3 ) 



Figure 2: (a) Lines of zero curvature R~ 1 (9s,k) = 0 (eqs. (3) and [10]) in the case of three- and 
two-dimensional drops (solid and dashed lines, respectively). ( 6 ) Parametric plot of R~ 1 (9B,k) 
(eqs. (4) [11] and [9]) for fixed three different values of 9b (the style of lines the same as in (a)). 
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In analogy to eq. [11] we also have 



(4) 


Solving eqs. (4) and [9] with respect to R we obtain R 1 = R l (9 


R 1 (9b, k), which is plotted for 


several values of 9b in Fig. 2b. 


III. SPHERICAL CAP. 

For the graphical definition of a spherical cap with an opening angle 9, see Figure 3. 



Figure 3: A plane intersecting a spherical drop divides it into two spherical caps - one (marked in 


grey) with an opening angle 9 and the other with an opening angle n — 6. 

IV. CALCULATION OF THE STABILITY DIAGRAMS FOR MULTIPLE JANUS 
DROPS. 

In the following we determine the shape of the liquid lenses forming a multiple Janus 
drop, given the ratio of volumes k = Vb/Va and the contact angles 9a and 9 B under the 
assumption that all the segments of the same phase are identical, i.e., that they form a 
periodic chain (section of the shortest repeating unit is sketched in Figure 4). 

First we determine the auxiliary angles [5a and /3b as functions of the radii Ra and Rb 
and the contact angles. The following relation holds: 


/5 a + [5 b + 9 a + 9 b — 2tt. 


(5) 


Additionally, from the law of sines for the triangle O 1 O 2 P (Fig. 4) we have 


Ra sin [5a = Rb sin /3b, 


( 6 ) 


which together with eq. (5) after applying basic trigonometric identities yields 

Rj sin(6U + 9 b ) 


(7) 


R L A Rj cos (9a + 9 b ) 
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Figure 4: Sketch of the section of the repeating unit of the multiple drop with graphical definition 
of the basic parameters. 
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Figure 5: The morphologic transition lines of a multiple drop for 0 = 1.15n calculated numerically 
using Maple (solid lines) and Surface Evolver (dotted lines). 

with i = A, B. The auxiliary angle (3 fulfills (see Figure 4) 

/3 = 9 b + ~ 7T- (8) 

The heights h* of the segments can be expressed as 


Jia/ 2 = Ra cos /3a~ R( 1 — cos /?), 

( 9 ) 

hs/ 2 = Rb cos /3b + R(l — cos/?), 

( 10 ) 
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where R can be either positive or negative. The condition of the balance of Laplace pressures 
yields the same result as for double Janus drops, i.e., (repeating eq. [8]) 

(sin 0 B sin6U\ -1 . la Q , f . 

R= [—R~ a - sm{e A + e B ). (11) 

Now fa, (3 and R are all functions of 9 A ,9 B , Ra, Rb- In order to change the independent 
variables from the unknown radii Ri to the given volumes V) one has to solve numerically 
the set of equations 

3Vk/(47r) = R 3 a ( 1 - 2/ 0 (ftO) - 2R 3 f 0 ((3), (12) 

3Vb/(47t) = i?|(l - 2/ 0 (ft,)) + 2R 3 fM (13) 

with respect to R A and R B , where fo(9) = (2 + cos 6) sin 2 (#/2). By choosing the length 
scale such that V A = 1 and V B = k one thus obtains fa = fa(9 A ,9 B ,k), fa = fa(9 A ,9 B ,k) 
and R = R(9 A ,9 B ,k). The conditions fa = 0 and fa = 7t/ 2 correspond to the limits of 
stability beyond which the segments must either deform or coalesce. For example, for a 
fixed value of 9 =: 0 A + 0 B the transition lines can be calculated in the (0 B ,k) plane by 
solving numerically the set of equations fa(9 — 0 B ,0 B ,k) = 0 and fa{9 — 0 B ,0 B ,k) = 7t/ 2 
for i = A,B. This yields four transition lines: k = fa hA \9 B ), k = fa hB \9 B ), k = k^ A \9 B ) 
and k = k^ 3 \9 B ). Depending on the value of 9 E [vr, 27 t] the lines and k^ j \ i ^ j, 

can either (i) intersect with each other (for 9 E \k,9*], where 9* is a certain angle which 

can be calculated numerically; we estimate 9* ~ 1.257T, see Fig. 4b in the main text) or (ii) 
intersect with the line 9j = n (for 9 E [9*, 2n]). In the hrst case (i) the continuation of the 
lines k^9 beyond the intersection points corresponds to the transition from the metastable 
to the unstable region. However, this continuation cannot be calculated within our model 
with spherical interfaces. The deformed shapes of the segments in this case can be calculated 
numerically by finite-element minimization of the interfacial energy using, e.g., the software 
Surface Evolver (see the dotted line in Fig. 5). 

We note that in Fig. 5 the lines /c^s) and fa 133 ' 1 meet at the single point (9 B = 7t/2, k = 0). 
This property can be actually proved for arbitrary 9. Assume that we approach the line k^ 3 ^ 
from the ’stable’ side, i.e., j3 B /* 7t/ 2 and simultaneously we take the limit k — V B —► 0 (we 
set V A = 1). Then from eq. (13) in which the left hand side vanishes and from the fact that 
/o(/?b /tt/2) 1/2 we have R 3 fo(/3 ) —> 0, which together with /o(/3) ~ (3 4 implies (3 —> 0. 
Next, we use another geometric relation: sin (3 = \R\~ 1 R B sm/3 B (see Fig. 4), from which it 
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Figure 6: (a) Top-view of a schematic representation of the nricrofluidic device with two nested 
T-junctions used for generation of double and multiple drops. Pressure sources (►), valves (XI) and 
inlets/outlets (o) to the chip (with the chip marked by dotted line) are indicated, (b) Side-view of 
the chip made of three layers of polycarbonate. The use of the middle layer enables positioning of 
the inlet to the chamber at the mid-height of the chamber, (c) Side-view of the chip made of two 
layers two layers of polycarbonate. In this case the chamber has the same depth as the channels. 

follows that /3 ~ \R\~ 1 Rb, which implies \R\~ 1 Rb —> 0. Now, we need to check if in this limit 
Rb remains finite. From eqs. ( 5 ) and (6) we obtain Rb = Ra sin(#i), where 6\ = 3n/2 — 9. 
Next from eq. (12) with Va = 1 we finally get R 3 B = 3sin 3 (# 1 )/[47r(l — 2f 0 (6 1 ))} 1 i.e., R B 
remains finite. This implies \R\~ X —> 0. Similarly, eq. (10) reduces to hs = R/3 2 ~ I R\~ 1 R 2 b 
which implies he — 1 > 0. This altogether shows that the lines k^ B \ k and the line R~ l = 0 
approach each other in the limit k — > 0. Finally, from eq. (8) it follows that in this limit we 
also have 6b —> tt/2. Analogous reasoning applies to the case k —> oo. 



V. MICROFLUIDIC DEVICES 

Schemes of the nricrofluidic devices are shown in Fig. 6. 
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